function out = CTBCS(x, y, z, r, b, r1, b1, a, b2, c)
%% CTBCS 系统迭代函数
    F = QuantumLogistic(x, y, z, r, b);
    G = SIMM(x, y, z, a, b2, c);

    arg_x = pi * (r1 * F(1) + (1 - r1) * G(1)) - b1;
    arg_y = pi * (r1 * F(2) + (1 - r1) * G(2)) - b1;
    arg_z = pi * (r1 * F(3) + (1 - r1) * G(3)) - b1;

    x_new = cos(arg_x);
    y_new = cos(arg_y);
    z_new = cos(arg_z);

    out = [x_new; y_new; z_new];
end